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ABSTRACT 

We have tested the apphcation to Sloan Digital Sky Survey data of the software package MATCH, 
which fits color-magnitude diagrams (CMDs) to estimate stellar population parameters and distances. 
These tests on a set of six globular clusters show that these techniques recover their known properties. 
New ways of using the CMD-fitting software enable us to deal with an extended distribution of stars 
along the line-of-sight, to constrain the overall properties of sparsely populated objects, and to detect 
the presence of stellar overdensities in wide-area surveys. We then also apply MATCH to CMDs for 
twelve recently discovered Milky Way satellites to derive in a uniform fashion their distances, ages 
and metallicities. While the majority of them appear consistent with a single stellar population, CVn 
I, UMa II, and Leo T exhibit (from SDSS data alone) a more complex history with multiple epochs 
of star formation. 

Subject headings: methods: numerical — methods: statistical — surveys — galaxies: dwarf — galaxies: 
stellar content 



1. INTRODUCTION 

Studying the properties of resolved stellar populations 
through color-magnitude diagrams (CMDs) can provide 
important constraints on stellar evolution and on the 
formation and evolution of stellar systems. For example, 
the apparent magnitude of the horizontal branch (HB) 
stars in a population, or the magnitude and color of 
the main sequence turn-off (MSTO), constrain or allow 
determination of the distance to and age and metal- 
licity of the population. The most detailed and most 
objectively quantifiable information can b e obtained by 
a numeric al CMD fi tting analys is (see e.g . Gallart et al.l 
1996; Tolstov fc Sa ha 1996; Ap aricio. Gall art fc Bertellil 



iDolphinI 119 97: Hohzm an et alj il999l : 



iHernandez, Gi lmore fc Valls-GabaudI 



Harris fc Z aritskv 2001). 

Over the last seven years, t he Sloan Digital Sky Sur- 
fSDSS. l\bAeral][2000l : lAdelman-McCarthv et all 
a large wide-field photometric 



vey 

|2007( ). has performed 



and spectroscopic survey in the northern Galactic cap. 
Because SDSS, unlike previous all-sky surveys, is a 
CCD-based survey, the photometric depth, accuracy, 
and homogeneity achieved is unprecedented over such 
a large area. SDSS data have proven to be of great 
value for the study of the Galactic stellar halo and 
nearby stellar systems with resolved popu l ations (e.g. 
Newberg ct al."2002; 'Belok urov et al.ll2006U lJuric et al.l 
2005; BeU et al. 2007). They have also led to the dis- 
covery of a host of p reviously unknown, faint Local 



Group members 
2006airbl: [B clokur 



( Willman et all l2005airbl: IZucker et all 
et al. I l2006d . I2007t llrwin et al.ll2007t 



Walsh et al.ii2007ll . 
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The study of dwarf galaxies and the Milky Way halo 
can provide important clues to the formation and as- 
sembly of the Local Group of galaxies and its individual 
constituents, as well as test cosmological theories. In 
the cold dark matter (CDM) paradigm, an abundance 
of subhaloes should be present around Milky Way type 
galaxies. The observed number of satellites is, however, 
one or two order s of magnitudes lower than what CDM 
theory predicts (|Klvpin et al.lll999t iMoore et al.|[l999l) . 
The recently discovered satellites partially alleviate this 
discrepancy, but if the CDM predictions are accurate, 
many more new satellite galaxies are to be discovered. 
Another prediction of CDM theories is that galaxies are 
formed hierar chically, i.e. built up by accreting smaller 
syste ms (e.g. ISearle fc ZinnI Il978t [Bullock & Johnsto^ 
l2005f) . Among the population of satellite galaxies sur- 
rounding the MW, some should have been formed within 
the MW halo while others should have a different origin. 
Studying the stellar populations, masses and abundance 
patterns of these satellites can provide important con- 
straints for galaxy formation and CDM theory. 

Our goals for this paper are two- fol d. First, we 
adapt the CMD-fitting software MATCH (jDolphinll200ll) 
for use on SDSS data, makin g use o f the theoretical 
isochrones provided by Girard i et al.l ([2002), and de- 
scribe several methods of applying this software to the 
data. A set of globular clusters taken from the SDSS 
data is used to test both the algorithmic approach and 
the theoretical isochrones. Second, we apply MATCH to 
the SDSS data for the new Milky Way satellites (that 
were recently discovered using the same SDSS data), to 
provide a uniform analysis of these objects and constrain 
their population parameters. 

In Section [2| we describe the data and the different 
CMD-fitting methods. The tests on globular clusters are 
done in Section [3| and in Section [3| we analyze the prop- 
erties of the new Galactic satellites. We conclude and 
summarize our results in Section [5] 



2 



de Jong et al. 



2. DATA AND METHODS 

2.1. Sloan Digital Sky Survey data 

The SDSS is a large photometric and spectroscopic 
survey, primarily covering th e North Galactic Cap. As 
of th e Data Release 5 (DR5, lAdelman-McCarthv et all 
[2001 . SDSS photometry covers almost a q uarter of the 
sky i n five passbands (u, q, r, i, and z; iGunn et al.l 
119981 120061 : iHogg et aLlboOlh . For the work presented 
in this paper we use photometry in the two most sen- 
sitive bands, namely g and r. A fully automated 
data reduction pipeline produces accurate astrometric 
and photometric measure ments for all d etected objects 
(iLupton. Gunn. fc Szalavl [19 99: Stough ton et al.l l2002l : 
Smith et al.l 120021 : iPier et all i2003^ Jvezic et al.l l2004 



Tucker et al.ll2006[ ). Photometric acc uracy is on the or- 
der of 2% down to g ~22.5 and r ~22 (jivezic et al.ll2004t 
ISesar et al. 2007.) . and comparison with deeper surveys 
has shown that the data are complete to r ~22. We use 
stellar data extracted from the SDSS archive with arti- 
facts remaved_aiid we_apply the extinction corrections 
from ISchlegel et all (|1998[ ) provided by DR5. 

In order for the CMD fitting with MATCH to work, a 
model describing the completeness and the photometric 
accuracy as a function of magnitude is needed. Com- 
pleteness has been checked for r-band'* by comparing de- 
tection of stars and galaxies by SDSS with the deeper 
COMBO-17 survey^. For stellar sources the complete- 
ness is 100% to r=22 and then drops to 0% at r=23.5. 
Since the stellar classification of COMBO-17 starts to be- 
come incomplete at approximately one magnitude fainter 
than in SDSS, this completeness is not extremely accu- 
rate, but amply sufficient for our purposes. The photo- 
metric errors given by the automated reduction pipelin e 
have been found to be very accurate (jivezic et al.ll2004f ). 
Comparison of the photometric accuracy in g and r as a 
function of apparent magnitude shows similar behavior, 
but in g the rise in the uncertainties occurs some 0.5 mag- 
nitudes fainter. For the g-band completeness we assume 
similar behavior as in r, but also shifted half a magni- 
tude fainter. Using this information we create artificial 
star test files that are given as input to our CMD fitting 
code. From this, the code creates the photometric accu- 
racy and completeness model with which the theoretical 
population models are then convolved. Figure [T] shows 
the completeness and average photometric errors of our 
artificial star files for the g- and r-band separately. Even 
though differences in, for example, seeing cause differ- 
ences in photometric accuracy and completeness between 
different parts of the sky, we use this same average model 
in the remainder of this paper. 

2.2. Color-Magnitude Diagram fitting 

Several software packages have been devel- 
oped for quantitative and al gorithmic CMD (or 
Hess diagram ) fitting fe.g . iGallart et all |1996l; 
Tolstoy fc S"ahal 119961 : Aparici o. Gallart & BcrtcUi 
1997t iDolphinI Il997t lloltzman et al l Il999i : lOlseni 
19991 iHernandez. Gilmorc & Vails- Gabaudl 1200 



Harris fc Zaritskvll200lD . The strength of these methods 
lies in the fact that all photometric information from 
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Fig. 1. — Completeness and photometric errors of the artificial 
star files adopted for the SDSS DR5 data and used for creating 
realistic stellar population models. Upper panel: g-band complete- 
ness fraction (solid line) and average photometric error (dashed 
line) as function of g-band magnitude. Lower panel: the same for 
r-band. 



imaging in two bands can be used simultaneously in a 
robust and numerically well-defined way. The software 
package we use is MATCH, which was written and 
successfully used for the CMD analysis of globular clus- 
ters and dwarf galaxies (|Dolphinl I2001D . The software 
works by converting an observed CMD into a Hess 
diagram (a two-dimensional histogram of stellar de nsity 
as a function of color and magnitude; iHessI 11923 ) and 
comparing it with synthetic Hess diagrams of model 
populations. To create these synthetic Hess diagrams, 
theoretical isochrones are convolved with a model of 
the photometric accuracy and completeness. The use 
of Hess diagrams enables a pixel-by-pixel comparison 
of observations to models and a maximum-likelihood 
technique is used to find the best-fitting single model, 
or best linear combination of models. To account for 
contamination by fore- and background stars, a "control 
field" CMD needs to be provided. MATCH will then 
use the control field Hess diagram as an extra 'model 
population' that can be scaled up or down. The resulting 
output will be the best-fitting linear combination of 
population models plus the control field. 

To construct star formation histories (SFHs) and 
age-metallicity relations (AMRs), a set of theoretical 
isochrones is needed that spans a sufficient range of 
ages, masses and metallicities. When MATCH was 
written, the only published set of isochrones with full 
age and metallicity coverage available was provided by 
iGirardi et al.l (|2002t l. Since we are specifically interested 
in the application to SDSS data, we need isochrones in 
the SDSS photo metric system. The o nly set available to 
date is one from IGirardi et all (|2004f l and has the same 
parameter space coverage. Some worries ex ist about th e 
accuracy of the isochrone calibration fe.g. IClemI 120061 ). 
which was limited to three systems available in the first 
data release (DRl) of SDSS. This unavoidably adds an 
extra source of uncertainty to all CMD fitting results. 
One of the main motivations for the various tests that 
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are described in this paper is to see whet her reasonable 
result s can be obtained using the current iGirardi et al.l 
(|2004l ) isochrones. 

2.3. Modes of CMD- fitting 

We use MATCH fo r slightly differ ent purposes than 
originally described bv lDolphinI (|2001h : the new "modes" 
are described briefly below. 

2.3. f. Star formation histories and metallicity 
evolution 

MATCH was originally developed to solve for the SFH 
and AMR of a stellar system for which all stars can be 
assumed to be at the same distance, such as globular clus- 
ters and dwarf galaxies. In this original mode of MATCH 
the distance and the foreground extinction to all stars 
are fixed a priori for each fit. From a set of age bins and 
metallicity bins, the combination that best fits the ob- 
served CMD is found, resulting in an estimate of the SFH 
and AMR. By doing the same fit with different values of 
the foreground extinction and distance, these two lat- 
ter parameters can also be determined by comparing the 
goodness-of-fit of the individual fits. In principle, param- 
eters like the steepness of the IMF and the binary fraction 
can be constrained in the same fashion, although in most 
cases the data quality and the degeneracies between dif- 
ferent parameters will prevent the determination of all 
of these parameters simultaneously. Analyzing the pa- 
rameter range of acceptable fits also gives robust error 
estimates for the SFH and AMR; this approach we refer 
to as MATCH'S "standard mode" . 

2.3.2. Distance solutions 

To apply CMD fitting techniques to nearby and ex- 
tended stellar syste ms (for e xample substructure within 
the Milky Way as in lde Jon g ct al. (2007)), the assump- 
tion that all stars are practically equidistant need not 
hold. To fit populations at differing distances simultane- 
ously we have implemented a new way in which to run 
MATCH that keeps the same number of free parameters, 
or parameter dimensions, to be explored. While in the 
standard mode the distance is fixed and age and metallic- 
ity are independent variables, here the distance becomes 
a free parameter, but age and metallicity are being lim- 
ited via a specified age metallicity relation (AMR). For 
any given AMR there is only one "population param- 
eter", e.g. age with a certain metallicity linked to each 
age bin. Of course, different age-metallicity relations can 
be explored, potentially covering a wide range of pos- 
sibilities. It seems sensible to choose an AMR that is 
appropriate for the studied stellar population on astro- 
physical grounds, such as the expected or known metal- 
licity evolution. For example, for the study of structures 
in the Galactic thin disk one would consider a different 
AMR than for the study of a metal-poor dwarf galaxy. 
However, within this new mode of running MATCH the 
choice of AMR is completely flexible and a flat AMR (no 
dependence of metallicity on age) or even an AMR where 
metallicity decreases with time is possible. 

2.3.3. Single Component fitting 

When studying low-luminosity or low surface bright- 
ness contrast systems, such as some of the recently dis- 
covered dSphs, or systems at low Galactic latitudes, for 



example substructures in the Milky Way, the contami- 
nation of the CMDs of these systems by fore- and back- 
ground stars will be very high. If the ratio of 'member' 
stars to field stars becomes low, the S/N of the CMD fea- 
tures to be fit consequently becomes very low. In such 
cases, all free parameters used to constrain distance, age, 
and metallicity in the MATCH mode of Sections 12.3.11 
and l2.3.2l can no longer be constrained by the data alone. 
To test a more restricted set of hypotheses in the case of 
low S /N systems we have developped a single component 
(SC) fitting mode for MATCH. In this mode, template 
stellar population models are created that consist of a 
single component with given, fixed distance, foreground 
extinction, age range and metallicity range^. Such a SC 
template is fit to the data together with a control field. 
Since the template consists of a single component, the 
only free parameters during the fit are the absolute scal- 
ing of (i.e. the number of stars in) the template and the 
control field. For each individual SC (i.e. for each combi- 
nation of distance, foreground extinction, age range and 
metallicity range), a value of the goodness-of-fit is ob- 
tained. Comparing the goodness-of-fit values for all the 
SC's we can then find the best-fitting templates and thus 
constrain the overall properties of the system. 

The quality of each fit is expressed by the goodness- 
of-fit parameter Q, determined using the max imum- 
likelihood technique described in iDolphinI ()2001l ) . This 
statistic is meant to measure how likely it is that the 
observed distribution of stars would result from a ran- 
dom drawing from the fitted model. In the SC fitting 
scheme we know that our fits will not be perfect since an 
SC model will not be a perfect representation of a real 
stellar population, but we are primarily interested in a 
likelihood-ratio test for a range of SC models. The ex- 
pected variance in the parameter Q for a certain CMD 
can be calculated from Poisson statistics and also ob- 
tained from Monte Carlo simulations using random draw- 
ings from the models. All fits with a value of Q less 
than f cr away from that of the best-fit SC model can 
be considered statistically acceptable, under the assump- 
tion that the best fit is a reasonable representation of the 
data. Taking all these fits together, error estimates on 
the fit parameters can be determined from the encoun- 
tered spread of the parameter values. 

This approach of using relative values of Q should only 
be used in cases where the "best fit" is a reasonable ap- 
proximation of the data. In cases where there are several 
very different prominent stellar populations this will not 
be the case and results will break down. Studying differ- 
ent parts of the CMD separat ely can help to circum vent 
such problems, as is done in Ide Jong et al.l ([2007t) and 
also in Section 14.1.11 of this paper for the dwarf galaxy 
Leo T. 

2.3.4. Over-density detection 

The same CMD fitting techniques can also be used 
to detect or verify, faint stellar over-densities in large 
datasets like that of the SDSS survey. In case of an 
overdensity at a specific distance, the stars will not only 
cluster around a certain spatial position, but also along 

^ While we use 'single components' it is important to note that 
these are not single stellar populations (SPP) with only a single 
age, distance, and chemical composition 



TABLE 1 

Artificial CMD properties 



CMD 


Ages (Gyr) 


[Fe/H] (dcx) 


m-M (mag) 


1 


10-13 


-1.5 


14.5 


2 


10-13 


-1.7 


14.5 




2-2.5 


-1.3 


14.5 


3 


10-13 


-1.5 


20.5 


4 


10-13 


-1.7 


20.5 




2-2.5 


-1.3 


20.5 



a certain locus in the CMD. The method we have devel- 
opped works by fitting the CMD of a small target patch 
of sky both with a control CMD obtained from a large 
area around the target patch and with the same control 
field plus a simple stellar population model. If no over- 
density is present, the fit with the control CMD should be 
good and no significant improvement in fit quality will 
result by including an extra stellar population model. 
However, if there is an over-density present and the pa- 
rameters of the population model, such as distance, age 
and metallicity, are close enough to those of the over- 
density, the fit including the model will result in a signif- 
icant improvement in the fit quality. By doing this for a 
range of models the presence of a stellar over-density can 
be assessed by looking at the difference in the fit quality 
parameter, AQ, between the control CMD fit and the 
control CMD plus model fit. Moreover, by comparing 
the values of AQ for different models, a rough estimate 
of the properties of the overdensity can be obtained at 
the same time. 

An important prerequisite for this method is of course 
that the over-density should contribute a larger fraction 
of the total number of stars to the target CMD than 
to the control CMD, and ideally should be contained 
completely within the target patch. This means that 
over-densities will be detected with the highest S/N if 
their size matches the size of the target patch. 

2.4. Fits to artificial data 

Before turning to real data it is useful to get a feel for 
the degree of precision that we can expect CMD fits to 
achieve. To do this, we create a few realistic artificial 
CMDs and apply SFH fitting, SC fitting and distance 
fitting methods to these to try and recover their proper- 
ties. U sing the theoretical isochrones fro mlGirardi et al] 
()2004D . a Salpeter initial mass function ()Salpeteilll955f ). 
and our model for the photometric errors and complete- 
ness of the SDSS data, we generate artificial stellar pop- 
ulations, to which field stars taken from the SDSS data 
are added. In total we create four CMDs, shown in Fig- 
ure m The number of stars in each CMD is chosen to be 
the same, to make for a fair comparison. CMDs 1 and 3 
contain a simple, single burst stellar population, located 
at different distances, so that either the MS and MSTO 
or the RGB and HB are prominent features. Artificial 
CMDs 2 and 4 contain two different age and metallic- 
ity populations, again with the CMDs corresponding to 
different distances. The detailed parameters of the arti- 
ficial populations are listed in Table [1] From the SDSS 
archives all stellar sources in a 5°x5° patch of sky close 
to the globular cluster M13 were extracted and a small 
fraction used to create the field star contamination in 
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Fig. 2. — Artificial color-magnitude diagrams used to test our 
CMD-fitting methods. Stars belonging to the artificial stellar pop- 
ulations listed in Table [T] are plotted thicker than the field stars 
taken from SDSS data. 

these CMDs. The remaining sources will be used as a 
"control field" population when fitting these background 
contaminated CMDs. 

We recover the SFHs of the artificial CMDs by running 
MATCH with nine age bins and twelve metallicity bins. 
The age bins are defined in log t because the isochrones 
are spaced more evenly in log t than in t, with the oldest 
running from 10.0 to 10.2 (10 Gyr to 15.8 Gyr) and the 
youngest from 8.6 to 8.8 (400 Myr to 630 Myr). The 
metallicity range of [Fe/H]= —2.4 to 0.0 is covered by 
twelve bins of 0.2 dex width. Thus, the artificial CMDs 
are fit with 108 model CMDs plus the "control field" 
stars, for a total of 109 partial CMDs. The best fitting 
linear combination of model CMDs will form the result- 
ing SFH and metallicity evolution. To get a sense of the 
uncertainties on the SFH, each CMD is fit thrice, each 
time with a slightly different assumed distance. The dis- 
tance moduli used are 14.3, 14.4 and 14.5 mag. The 
resulting SFHs and metallicities are presented in Figure 
[3]as the relative star formation rate (SFR) and as [Fe/H] 
as function of the age of the stars. 

For the SC fits models are created for a grid of pa- 
rameters with a distance modulus range of 2.0 with 0.1 
magnitude intervals centered on the the actual distance 
modulus (Table [T]). Ages used run from 9.0 to 10.2 in 
log t with bin widths of 0.1 dex and metallicities go from 
[Fe/H]= —2.1 to 0.0 in steps of 0.1 dex and bin widths of 
0.2 dex. Each model plus the "control field" is fit to an 
artificial CMD, after which the resulting goodness-of-fit 
values are compared. The resulting contour plots show- 
ing the best-fitting parameter combinations are shown in 
Figure [H 

Using the new distance solution option of MATCH we 
solve for the SFH of the artificial CMDs again while fix- 
ing the metallicity to a single, uniform value. The same 
age bins are used as for the SFH fits described above, and 
a distance range of 3.0 is probed with 0.1 magnitude res- 
olution. Including the "control field" , 271 partial CMDs 
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Fig. 3. — Star formation histories for the artificial CMDs from the "standard mode" of CMD-fitting. In the columns on the left the SFR 
relative to the average SFR over the total age range is plotted, and in the columns on the right the SFR-weightcd metallicity for the bins 
with a significant SFR. Horizontal error bars indicate the width of the age bins. The vertical error bars on the SFR give the full range of 
SFRs found in the fits, while the vertical error bars on the metallicities give the standard deviation in the metallicities of the fits. Dashed 
lines indicate the input SFHs and metallicities for the artificial CMDs. 



are used in these fits. The fits are done at three different 
values of the metallicity: [Fe/H]= —1.6, —1.5, and —1.4 
to probe the uncertainties. Figure O shows the resulting 
SFHs in relative SFRs and distance determinations. 

The SFHs in Figures [3] and [5] show that the software 
succesfuUy recovers the single burst in CMDs 1 and 2, 
and the two distinct episodes of star formation in CMDs 

3 and 4. Some bleeding of star formation in adjacent bins 
will occur even when the fitted age bins exactly match 
the input SFH (J)olphin 2001^. Therefore it is to be 
expected that this also happens here, since the duration 
of the input star forming episodes do not match any age 
bins perfectly. It should be noted that the error bars 
in the SFHs represent the range of SFRs encountered in 
each age bin within the three fits on which the SFHs 
are based, which means that the errors between adjacent 
bins are correllated. For example, in the SFH of CMD 

4 in Figure [3] the two bins at ages of 2 to 4 Gyr are not 
independent; when the SFR in one of them is low, it 
is high in the other. That the age found for the young 
population in CMD 4 in Figure [5] is the natural result 
of the metallicities in the fit being lower than the actual 
metallicity of the younger population. In all cases, the 
average metallicities and distances found in Figures[3]and 
[5] agree well with the input values. 

Turning to the SC fits in Figure HI we see that for 
CMDs 1 and 3, the input properties (indicated by stars) 
are recovered very well. Because there is a degeneracy 
between distance, age and metallicity when only the MS 



and MSTO are seen, the properties of CMD 1 are, how- 
ever, constrained more poorly than for CMD 3. Since 
the SC fits use a single model population, the results 
for CMDs with multiple populations, such as CMDs 2 
and 4, cannot be expected to properly recover the prop- 
erties of all stars. Rather, the population representing 
the majority of stars in the CMD will have the strongest 
weight. This is clearly visible in the SC fitting results 
for CMD 4, that recover the properties of the older pop- 
ulation, that dominates the star counts in the RGB and 
HB. For CMD 2, however, the upper MS and MSTO of 
the younger population has such high S /N that this pop- 
ulation does influence the SC fits significantly and the 
favored age is relatively young ('^4 Gyr). This shows 
that for high S/N CMDs with a complex stellar popula- 
tion make-up, the SC fits is not well suited; they were, 
however, developed for application to low S/N CMDs. 

What distinguishes these artificial CMDs from the real 
data, is that the artificial stellar populations are based 
on the exact same isochrones, IMF, and photometric er- 
ror and completeness model as used by MATCH. Inac- 
curacies in any of these will degrade the fit results one 
obtains for real data compared to the results for the ar- 
tificial data. 

3. TESTS ON GLOBULAR CLUSTERS 

Globular clusters arc ideal test cases for the SDSS im- 
plementation of MATCH since in most of them their stars 
have been found to be of a single age and nearly single 



6 



de Jong et al. 











r 












BJZ 






BLS B.^f BL.fi i.a iru] 




-M DUO 



Id 






Fig. 4. — Contour plots of the single component fitting results for the artificial CMDs. From left to right, the columns show: distance vs. 
metallicity; distance vs. age; age vs. metallicity; the CMD with the best-fit single component model in contours. The grey-scale contours 
show the regions in parameter space that are within la, within 2a, within 3a and outside 3a of the best-fitting single component model, 
going from white to dark gray. Input values for the artificial stellar populations (see also Table [TJ are indicated by asterisks. 



metallicity. Here we will use six GCs present in the SDSS 
data to test how we ll MATCH, using the isochrones from 
iGirardi et aDlpOOl ). can recover their known population 
properties. The systems in question, M13, M15, M92, 
NGC 2419, NGC 6229, and Pal 14, cover a large range 
in distance, background contrast (i.e. ratio of member 
stars versus field stars) and have metallicities between 
[Fe/H]= -2.3 and -1.3 dex. For the following tests, all 
stars within a 5°x5° field centered on each GC were 
obtained from the SDSS catalog archive server (CAS). It 
should be noted that in the cases of M15, M92 and NGC 
2419 these fields were not completely covered by the DR5 
footprint. Flags were used to exclude objects with poor 
photometry^ and all stars were individually corrected for 
foreground extinction using the correction fact ors pro- 
vided by DR5 and based on the maps by Schle gel et al.l 
([1998). Any contamination of the samples of stars ob- 
tained by misclassified galaxies or other contaminants do 
not pose problems for our analysis because of the use of 
control fields that should contain the same populations 
of these contaminants. "Control field" CMDs were con- 
structed from all stars outside a radius of 0.5° from the 

See |http://cas.sdss.org/astro/en/help/docs/realquery.asp#flags| 



center of each GC. "Target fields" were constructed by 
selecting stars within a certain radius from the center. 
The radu used for M13, M15, M92, NGC 2419, NGC 
6229, and Pal 14 are 0.25°, 0.2°, 0.25°, 0.15°, 0.15°, and 
0.1°, respectively. These radii are similar to the tidal 
radii of the GCs but are not strictly related to them, 
since in all clusters - and especially in the densest (M13 
and M15) - the SDSS photometry in the central parts 
is missing due to the high degree of crowding. With the 
radii used, a significant number of field stars is present 
in each target CMD, so that MATCH has a good handle 
on the proper scaling of the control CMD. 

The SDSS CMDs of the GCs are shown in Figure M 
and their previously measured parameters (distance, age, 
[Fc/H]) are listed in Table [21 the mean reddening M^B — 
V) in magnitudes according to the Schlegel e t al.l ()1998l ) 
maps in the target area for each GC is also listed in the 
table. The GCs in the three top panels in Figure [6] have 
metallicities of around [Fe/H]= —1.5 dex, and the GCs 
in the lower panels are more metal-poor ([Fe/H]< —2.0). 
For the more nearby systems (M13, M15, M92 and NGC 
6229) the main CMD features that will drive the fits 
are the main sequence (MS) and main sequence turn off 
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Fig. 5. — Star formation histories and distances for the artificial CMDs from the "distance solution" mode. The columns on the left 
and the error bars are as in Figure |3] In the columns on the right now the SFR- weighted average distance is plotted for the age bins with 
significant star formation. Dashed lines indicate the actual input distance for the artificial CMDs. 



TABLE 2 

Globular cluster properties 



Cluster 


R© (kpc) 


m-M (mag) 


[Fe/H] (dex) 


Age (Gyr) 


< E{B - V) 


M 13 


7.7, 7.5") 


14.38, 14.43 


-1.54 


12°) 


0.019 


M 15 


10.3, 10.0''' 


15.37, 15.0'') 


-2.26 


13.2'') 


0.11 


M 92 


8.2 


14.6 


-2.28 


14c,d) 


0.025 


NGC 2419 


84.2, 98.6'=) 


19.63, 19.97'=) 


-2.12 


14'=) 


0.060 


NGC 6229 


30.4 


17.44 


-1.43, -1.1-f) 


~11.59) 


0.023 


Pal 14 


73.9 


19.34 


-1.52 


~10'^) 


0.035 



Note. — Values ar e from Harris (1996), ex cept fo r: a) Grundahl et al. (1998); b) 
IMcNamara et al.l (l200l): c'jlPoiit et al.,, (,19981): d) [Grundahl et al. 12000) : e) Harris et"an 
(119971 ): D IBorissova etaLl a999l) : e) IVandenbera I2OOOI) : h) ISaraiedn?n997l) 



(MSTO); for the more distant systems (NGC 2419 and 
Pal 14) only the red giant branch (RGB) and horizontal 
branch (HB) are observed. Pal 14 in particular represents 
the case of a very sparsely populated system with high 
contamination from field stars. 

In the following subsections, we discuss the results of 
applying MATCH to these GCs in the four different 
modes described in section 12.31 In fitting the CMDs 
of these objects we will use a color range of —0.5 < 
g — r < 1.5, and limit the magnitude range to g < 22.5 
and r < 22.0; the color cut includes the MSTO and the 
tip of the RGB, but excludes a large fraction of the nu- 
merous low-mass foreground stars. For all the fits in this 
section we will assume a slope for the stellar initial mass 
function that is close to Salpeter; the binary fraction is 



assumed to be 0.1, in line with recent observational and 
theor etical evidence for low binary fractions in globulars 
(e.g.. lIvanova et al.ll2005t IZhao &: Bailvnll2005f) . 

3.1. Star formation histories 

To measure the SFH and AMR of the clusters we run 
MATCH fits with nine age bins and twelve metallicity 
bins. Because the isochrones are spaced more evenly in 
log t than in t, the age bins are defined to have equal 
width in log t, with the oldest running from 10.0 to 10.2 
(10 Gyr to 15.8 Gyr) and the youngest from 8.6 to 8.8 
(400 Myr to 630 Myr). The metallicity bins cover the 
range from [Fe/H]= —2.4 to 0.0 in twelve bins of 0.2 dex 
width. Together with the control CMD this yields 109 
model CMDs for each GC. The fitting software then finds 
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Fig. 6. — Color-magnitude diagrams (CMDs) from SDSS DR5 data for the six globular clusters with independently determined population 
parameters, used to test the CMD-fitting software. 



the linear combination of these model CMDs that best 
fits the observed CMD. Each object is fit several times 
with slightly different fixed distances to determine the 
impact of distance errors on the recovered solutions. In 
the cases of M15 and NGC 2419, where values from the 
literature for the distance modulus are quite discrepant, 
we run fits for distance moduli of 15.0 to 15.4 (M15) and 
19.6 to 20.0 (NGC 2419) in 0.1 magnitude steps. For the 
other clusters we run three fits each, one at the literature 
value and at 0.1 magnitude further and nearer. 

Figure [7] shows the results of these fits, with the frac- 
tional SFR as a function of time in the left column, and 
the metallicity as a function of time in the right column. 
Metallicities are only plotted for age bins with a signif- 
icant SFR. The plotted SFR for a given age bin is the 
average of the values given by the fits at the different dis- 
tances, with the error bars showing the complete range 
of SFR values found. The metallicity in a bin is the av- 
erage of the fits weighted by the SFR, and the error bars 
correspond to the standard deviation. 

Looking at the three intermediate metallicity GCs in 
the top half of Figure [7] we see that in all cases MATCH 
infers a purely old (>7 Gyr) population, as expected. For 
M13 and NGC 6229 practically all stars are found to be 
older than 10 Gyrs; some bleeding into the second old- 
est age b in is always likely to occur, cf. Section [2T4l and 
iDolphinl (j200l|. In the case of Pal 14 the SFR is found 
to be constant over the two oldest age bins, but this onl y 
confirms the age estimate of '^10 Gyr (jSaraiedinilll997t ). 
which implies that star formation would likely have been 
spread over these two bins. The intermediate metallici- 
ties of M13, NGC 6229 and Pal 14 are reproduced weU 
and our recovered values are within 0.1 to 0.2 dex of 
the literature values. For N GC 6229 our f its favor the 
slightly lower metallic ity from [ Harris! (1 19961) rather than 
the higher value from iBorissova et al.l 

Turning to the bottom half of Figure [7| with the three 
metal-poor clusters, we see that here also the results are 



consistent with a purely old population, although the 
SFHs of M15 and NGC 2419 show a non-zero SFR in 
a few younger age bins. In the case of M15 this is due 
to the blue horizontal branch (BHB) stars that are not 
properly reproduced by the theoretical isochrones, caus- 
ing the code to fit them as a younger turn-off. This 
demonstrates that CMD fitting results always need to be 
looked at in detail in order to understand their implica- 
tions. The age distribution in NGC 2419 is not very well 
constrained because no MSTO is present in the CMD. 
The metallicities of M15 and M92 are recovered well for 
the old stars. While the estimated metallicity of —1.9 
for NGC 2419 is on the high side, for these metal-poor 
GCs the recovered metallicities are within 0.2 dex of the 
literature values. 

3.2. Single component fitting 

As described in sect ion [2 . 3 . 31 above . in the Single Com- 
ponent (SC) fitting mode an individual model CMD is 
fit to the data together with the control CMD, rather 
than a whole set of models. The results therefore pro- 
vide information on the characteristic properties of the 
dominant stellar population in the CMD. We sample the 
3D parameter space of distance, age and metallicity by 
using one of nine age bins, one of 21 distance bins, and 
one of 23 metallicity bins. The age bins have a width 
of 0.1 log-years, going from 10.1-10.2 log-years (12.5 Gyr 
to 15.8 Gyr) to 9.7-9.8 log-years (5 Gyr to 6.3 Gyr) in 
steps of 0.05 log-years. This age range more than suffices 
for these GCs. We use a distance modulus range of two 
magnitudes centered on the literature value with a res- 
olution of 0.1 magnitude. Metallicity bins are centered 
at [Fe/H]= —2.2 to 0.0 with an inter-bin distance of 0.1 
dex and a width of 0.2 dex. 

After matching each of the 4347 SC options to each GC 
and determining the expected variance in the goodness- 
of-fit parameter Q, we take the best fit as our measure- 
ment and compare the other fits with respect to these 
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Fig. 7. — Star formation histories and age metallicity relations for the six globulars obtained by using MATCH in its "standard" mode 
(fitting for age and metallicity distributions at a given distance). Star formation rates are shown relative to the average star formation 
rate over the entire age range. Error bars along the age axis show the width of each age bin, while the error bars along the y-axis show 
the range of star formations in the age bin at the three different distances for which a fit was done and the uncertainty in the metallicity 
in each age bin. Metallicities are only shown for bins with significant star formation. The dashed lines in the right panels indicate the 
literature values for the metallicity from Table |2] 



best values. The resulting likelihood contour plots for 
metallicity, distance and age measurements are shown in 
Figure [SI Note that in these contour plots, the third pa- 
rameter is unconstrained, rather than fixed at its known 
value. If the third parameter would be kept fixed, the 
contours would shrink significantly. In general the agree- 
ment with the literature values (indicated with asterisks 
m Figure [H) is very good. The small offset for M13 is 
most likely due to inaccuracies in the isochrones, which 
illustrates the added uncertainty in CMD fitting results 
that depend on a single set of theoretical isochrones. Our 
distance modulus measurement from the SC fits for M15 
lies exactly in between the two literature values, while for 
NGC 2419 our results c l early favor the smaller distance 
value from [Harris et all ()1997f ). From the two metallic- 
ity determinations in the literature for NGC 6229, the fit 
results agree significantly better with the lower metallic- 



ity of [Fe/H]^ —1.4 dex. In some cases, it looks as if 
the peak of the likelihood lies outside of the parameter 
space probed, at lower metallicities or older ages. Unfor- 
tunately, these limits of the parameter space are dictated 
by the parameter range of the available isochrones. 

Also shown in the right column of Figure [5] are the ob- 
served CMDs of the GCs, with contours representing the 
model populations. In most cases we used the literature 
parameters for these models, except for M15, where we 
chose a distance modulus of 15.2, the mean of the two 
different available values. For NGC 6229 and NGC 2419 
we use the literature metallicities and distance moduli 
that best correspond to our results. The agreement of 
model and data is generally acceptable, although inaccu- 
racies, especially in the sub-giant branch, can be seen in 
M13 and M15. 

Important to note is the influence of the degeneracies 
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Fig. 8. — Comparison of single component population parameters from SDSS and MATCH versus their literature values for the six 
globulars. The columns show (from left to right): distance vs. metallicity; distance vs. age; age vs. mctallieity; the color-magnitude 
diagram with the best- fit single component model in contours. The grey-scale contours show the regions in parameter space that are within 
1(7, within 2<T, within 3a and outside 3cr of the best-fitting single component model, going from white to dark gray. The solid line contours 
are for the case in which one of the parameters is fixed to its best-fit value, i.e. age for the left-most, metallicity for the middle-left, and 
distance for the middle-right columns. Literature values for the parameters are indicated by asterisks. 
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Fig . 9. — Star formation histories and distances for the six globulars obtained by using MATCH in the distance solution mode (Section 
12. 3. 21 1. As in Figure [7] star formation rates are shown relative to the average star formation rate over the entire age range and horizontal 
error bars show the width of the age bins. Vertical error bars indicate the range of star formation rates encountered at the three different 
metallicities for which a fit was done, and the uncertainty of the distance. The dashed lines in the right panels indicate the distance values 
from the literature (see Table [J}- 



between age, distance, and metallicity. For the nearby 
systems, M13, M15 and M92, where the fit is driven by 
the MS, there is a clear distance- metallicity covariance; 
the left-most panels in Figure [8] show that the distance 
estimate increases with increasing metallicity. For the 
distant GC Pal 14, where the fit is driven by the RGB, 
the covariance is the other way around. Another impor- 
tant, yet not wholly surprising, point is that higher S/N 
in the CMD features results in much tighter contours, il- 
lustrated by the fact that M13, M15 and NGC 2419 give 
the best constrained results. NGC 2419 has the addi- 
tional advantage that the well-populated HB provides a 
very strong constraint on the distance, making this sys- 
tem the best-constrained of the lot. 

3.3. Distance solutions 

As mentioned in Section [2.3.21 we explore fits for dif- 
ferent distances and different ages, with the metallicity 



linked to the age. The result of each fit is a SFH with a 
(potentially different) distance for each age bin; differing 
distances for different sub-populations are not astrophys- 
ically relevant for globular clusters, but serve as a test 
case. For fitting the GCs it seems reasonable to use the 
same metallicity in all age bins. We use the same nine 
age bins as in Section 13.11 In distance space we use a 
distance modulus range of three magnitudes centered on 
the known distance of each system with a resolution of 
0.1 magnitudes, giving a total of 31 distance bins. For 
each GC we carry out three fits with different metallici- 
ties, centered on the literature values (see Tabled with 
offsets of 0.1 dex. For NGC 6229 we assume the lower 
metallicity value [Fe/H]= —1.43, which is favored by our 
SFH and single component fitting results. 

In Figure [9] we present the age-distance distributions 
that result from these fits. Compared to the fits in sec- 
tion l3. li the number of free parameters per fit is the same, 





Fig. 10.— match analysis of a composite CMD: aj CMD of the globular M 13. fij CMD for the globular NGC 6229. cj combined CMD 
of the globulars M 13 and NGC 6229. d) average result of the three distance solution fits to the CMD in panel a), showing the number of 
stars in each age-distance bin; the M 13 stars are fit by the oldest age bins around m-M=14.4 magnitudes, e) as in d) but for the CMD in 
panel b); the NGC 6229 stars are fit by the oldest age bins at m-M=17.2 and 17.4 magnitudes, f) as in d) and e) but for the fits to the 
combined CMD in panel c); both the ages and distances of the two distinct populations are correctly recovered. 



with now metallicity fixed instead of distance. Thus, 
in general similar quality SFHs should be expected, al- 
though not necessarily exactly the same. While in the 
"standard" SFH fits the star formation within an age 
bin can be spread over several metallicity bins, it can 
now be spread over several distance bins. These spread 
will translate diff'erently into the color-magnitude plane, 
and therefore can result in differences between the SFHs 
in Figures [7] and [9] 

As in Figure [7] all GCs are found to have a single peak 
in their SFR in the oldest age bin, in some cases with 
some bleeding into the adjacent bin. The extension of 
the SFH of Pal 14 to more recent times has large uncer- 
tainties and is not significant. It should be taken into ac- 
count that the exact age of this system is difficult to con- 
strain with SDSS data, since the MSTO is not present. 
Interesting to see is that M13 and M92 show almost no 
bleeding into the second oldest age bin here, while they 
did in Figure [71 for M15 and NGC 6229 the exact oppo- 
site is the case. Apparently, the former are fit better with 
some distance spread, while the latter benefit more from 
a metallicity spread. This does not necessarily reflect an 
actual distance or metallicity spread, but rather shows 
that inaccuracies in the isochrones are sometimes better 
accounted for one way and sometimes another way. 

The distances we find for the oldest age bin, and there- 
fore for the majority of the stars, show very good agree- 
ment with the literature values. Comparing Figures [7] 
and [5] shows that the SFHs obtained when fixing either 
the distance or the metallicity are very similar. 

The main motivation for enabling this distance fitting 
mode was, however, not to analyze single-distance sys- 
tems, but to be able to cope with cases where stars 
are located at diff erent distances (cf. Canis Major, 
Ide Jong et al.l[2007h . To demonstrate this apphcation, 
we combine the CMDs of M13 and NGC 6229, resulting 



in a CMD containing two populations with similar metal- 
licities and ages, but located at significantly different dis- 
tances. The CMDs of M13, NGC 6229 and the combined 
CMD are shown in panels a), b) and c) of Figure fTOl We 
run distance solution fits to these CMDs in the same 
way as before, with the same age bins, but now using a 
larger distance modulus range, from 12.0 to 19.0 magni- 
tudes with a resolution of 0.2 magnitudes. Three fits are 
run, for uniform metallicities of [Fe/H]= —1.4, —1.5, and 
— 1.6. We include also a control CMD in the fit which is 
the control CMD of M 13; since M 13 and NGC 6229 are 
only 10 degrees apart and at similar Galactic latitudes, 
this gives good results. The averaged best-fitting distri- 
butions of stars over the grid of age and distance bins are 
shown in panels d), e) and f) of Figure [TUl Both ages and 
distances of the two populations are recovered extremely 
well, with all stars in the oldest age bin (>10 Gyr), M13 
stars centered at a distance modulus of 14.4 and NGC 
6229 stars at a distance modulus of 17.4. Moreover, the 
result for the combined CMD is almost exactly the sum 
of the results for the individual CMDs, clearly illustrat- 
ing the power of MATCH to decompose more complex 
CMDs. 

3.4. Over-density detection 

The final application we describe and test here is aimed 
at algorithmically detecting stellar populations which are 
confined to a limited distance and area on the sky, i.e. 
that form "overdensities" . Of course, the GCs used in 
this section are very obvious stellar overdensities and 
should be detected at very high S/N. Following the ap- 
proach described in section 12.3.41 we define a grid of 
20'x20' boxes in each 5°x5° field that was obtained 
from the SDSS catalog server, with the central box cen- 
tered on the GC. The CMD of each box is fit with the 
control field CMD and with the control field plus a sim- 
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Fig. 11. — Stellar overdensity detection using MATCH: each panel shows the fit quality improvement, dQ, between fitting only a control 
field and fitting a control field plus stellar population model to the CMDs of 20" X 20" regions around the six globular clusters. The fields 
plotted are 5°x5°fields centered on the globulars with right ascension along the x-axis. 



pie population model. For each GC we use a model with 
the appropriate distance, age and metallicity. Figure [TT] 
shows the difference in the goodness-of-fit parameter Q 
between the two fits for areas around each GC. Note that 
for M15 and M92, the 5°x5° field is not completely ob- 
served by SDSS. Clearly, all GCs produce high peaks in 
AQ. Since AQ is roughly proportional to the number 
of stars attributed to the population model in the model 
fit, the objects with the largest numbers of stars produce 
the strongest detections. Only in the plot for Pal 14 does 
the vertical scale allow the random variations in AQ to 
be visible in the area around the object. The S/N of 
these detections should be measured with respect to this 
random noise. Using the 216 boxes around the central 
9 boxes we obtain an estimate of the rms in AQ from 



random variations in the CMDs. Negative values of AQ 
are suppressed because a negative number of stars in a 
population is unphysical, which has to be taken into ac- 
count when calculating the rms. Even Pal 14, a rather 
sparsely populated and distant cluster, is detected with 
high S/N, as AQpaiu/rmsAQ — 56. 



4. APPLICATION TO NEW MILKY WAY 
SATELLITES 

In summary, Section [3] has illustrated the robustness 
and precision with which MATCH can interpret high- 
quality CMDs in a variety of circumstances. Over the 
past few years the importance of deep wide-field sur- 
veys for the study of the stellar components of the 
Milky Way gala xy has been de monstra ted by 2MASS 
(e.g.,,Oiha,2qo iHlbata et al.ll2002: Majewski et all l2003l: 
iBrown et al.l l2004: Martin et alJ I2004D and the SPSS. 
Specifically, the SDSS has enab led a detailed picture of 
the Galactic st ellar halo (e.g., iBelokurov et al.l [2006bl : 
iBell et al.l l2007l and led to the discovery of at least four- 
teen previously unknown d warf galaxies an d globulars: 
Ursa Maior I (UMal, Will man et al.ll2005a[ ). Willman 1 
CWill. IWillman et al.i ,2005bl). Canes V enatici I (CVnl 
Zucker et al.l l2006af) . Bootes I f BooI. IBelokurov et alj 
2006d), Ursa Major II (UMall, iZucker et al. l l2006bD . 
Coma Berenices, Canes Venatici II, H ercules, Leo IV, 
Segu e 1 (Com, C Vnll, Her, LeoIV, Segl, IBelokurov et al.l 
20071) ■ Leo T lllrwin et al.l 120071 ). Bootes II (BooII, 
Walsh et al1l2007t l. and two very low-luminosity globu - 
lars Koposov 1 and Koposov 2 (|Koposov et al.1 [2007aD . 
While these discoveries are important for our under- 
stand ing of the so-called "missing satellit e problem" 
(e.g.. iKlvpin et al.lll999l: iMoore et all 119990 . these new 
systems might also shed new light on the mechanisms of 
star formation in low-mass halos with small baryon con- 
tent. Why these systems were only discovered with SDSS 
data can be explained by their low surface brightness, 
as shown by iKoposov et al.l (|2007bD . Apart from the 
promise that more systems of even lower surface bright- 
ness might be waiting to be discovered by future genera- 
tions of deep imaging surveys, this fact raises intriguing 
questions. Do these systems show that star formation 
can occur under conditions generally considered to be 
incapable of supporting such activity, or have these sys- 
tems evolved to their current state through, for example, 
tidal interactions with the Galaxy? A thorough and sys- 
tematic study of the properties of these objects is needed 
to understand their formation and evolution and the pos- 
sible cosmological implications. 

Several groups have been working on photometric and 
spectroscopic follow-up of the new MW companions in 
order to constrain their properties. A compilation of 
these results is given in Table |3l Although individual 
measurements of properties might be less precise in some 
cases, the advantage of studying these objects using the 
SDSS photometry is that this highly uniform dataset al- 
lows a systematic approach to the complete set of objects 
and fair comparisons between them. In the following, we 
apply our single component (SC) fitting techniques to 
constrain the overall distances, metallicities and ages of 
the objects listed in Table |31 We also constrain the SFH 
of a subset of objects, since the majority contain too few 
stars to make this a worthwhile exercise. The objects we 
use are Boo, CVn I, and UMa I because of their relatively 
well-populated CMDs, and Leo T because of its clearly 
complex star formation history. For all the fits in this 
section we assume a Salpeter-like initial mass function 
and a binary fraction of 0.5, similar to the value for d isk 
stars (jDuauennov fc MavoHri99lUKroupa et al.lll993f ). 



4.1. Fitting methods and results 
4.1.1. Single component fits 

As for the GCs we obtained data in 5°x5° fields 
around each object from the SDSS CAS as outlined in 
section [3] For the target CMDs we used the stars within 
one half-light radius (r/i, Table[3|) from the center (with a 
minimum of 3'), with all stars outside a radius of 60' used 
for the control CMDs. Figure [T2l shows the target CMDs 
for the twelve objects. In all cases except Leo T a color 
range of —0.5 < g — r < 1.5 is used for the CMD fits. The 
population parameter space we sample consists of 12 age 
bins, 24 metallicity bins and 31 distance bins. The age 
bins have a 0.1 log- year width and separations of 0.1 log- 
year, with the oldest running from 10.1-10.2 log-years 
(12.5 Gyr to 15.8 Gyr) and the youngest from 9.0-9.1 (1 
Gyr to 1.26 Gyr). Metallicities run from [Fe/H]=-2.3 
to 0.0 in steps of 0.1 dex, and each metallicity bin has a 
width of 0.2 dex. The range in distance modulus probed 
is 3 magnitudes wide, with a resolution of 0.1 magni- 
tudes, and is centered on the literature value for each 
object. 

For Leo T we use a slightly different approach. The 
group of blue stars in the CMD is too bright and does 
not have the right morphology to be interpreted as HB 
stars belonging to the same population as the RGB stars. 
This means there are two populations present in two dif- 
ferent parts of the CMD, with the RGB corresponding to 
the old population and the stars at g — r <0.3 to a much 
younger population. As demonstrated in Section [2.41 in 
such a case the SC fits will give very poor results. How- 
ever, since the two populations in Leo T are only present 
on one side of the g — r ~ 0.3 line, they can be easily 
separated. Thus, to determine the properties of the two 
populations separately we divide the Leo T CMD into a 
red (-1.0< g-r <0.3) and a blue (0.3< g-r <1.5) part 
and fit these separately. The red part of the CMD is fit 
with the same parameter space as the other objects. For 
the blue part the same metallicities and distance moduli 
are probed, but a younger set of age bins is used. Here 
we use nine bins, the four oldest of which have a width 
of 0.1 log-years centered on 9.05, 9.15, 9.25 and 9.35; the 
five youngest have a width of 0.2 log- years centered on 
8.1, 8.3, 8.5, 8.7 and 8.9 log-years. 

In Figures [T3l and [141 the SC fitting results for the new 
MW satellites, except Leo T, are presented. The right- 
most panels show the CMD of each object again, with the 
best-fit SC model overplotted as contours. Monte Carlo 
runs were used to determine the values of a. Except for 
Boo I and CVn I, solid contours are shown for the case 
of a fixed metallicity or age, where such a measurement 
is available in the literature. Figure [15] shows the same 
results for the red and blue part of the Leo T CMD. The 
fact that low S/N hinders a robust determination of the 
population properties provides a serious problem for the 
majority of the targets. Although for most objects a fi- 
nite region of parameter space lies within 1 a from the 
best fit (the white area in the contour plots), reasonable 
measurements of any individual parameter can only be 
obtained for Boo I, CVn I, CVn II, Com and Her. Espe- 
cially for Boo II, Leo IV and Leo T, parameters are very 
poorly constrained, and Q space shows a complicated 
morphology with multiple minima. 
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Fig. 12. — Color-magnitude diagrams from SDSS data for twelve recently discovered Milky Way satellites. The dashed lin e in th e panel 
for Leo T indicates the separation between the red and blue parts used in the single component fitting analysis (see Section l4.1.1ll . 



TABLE 3 
New MW satellite properties 



Object 


Rq (kpc) 


m-M (mag) 


[Fe/H] (dex) 




Mv (mag) 


Bootes I 


60±6 


18.9±0.2, 19.1±0.1'"\ ig.OztO.l') 


~-2.5'=', ~-2.1'') 


12'.8±0'.7 


-5.8±0.5 


Bootes II 


60±10 


18.9±0.5 




4'±2' 


-3±1 


Canes Venatici I 


220±20 


21.75±0.20 


-2.5 - -1.5=', -2.0'*'. -2.09^' 


8'.5±0'.5 


-7.9±0.5 


Canes Venatici II 


150±15 


20.9±0.2 


-2.3-f) 


~3'.2 


-4.8±0.6 


Coma Berenices 


44±4 


18.2±0.2 


~-2, -2.00-'') 


~5'.5 


-3.7±0.6 


Hercules 


140±13 


20.7±0.2 


-2.27^^) 


~8'.2 


-6.0±0.6 


Leo IV 


160±15 


21.0±0.2 


-2.3/) 


~3'.4 


-5.1±0.6 


Leo T 


~420 


--23.1 


~-1.6, -2.3-') 


~1'.4 


~-7.1 


Segue 1 


23±2 


16.8±0.2 




~4'.5 


-3.0±0.6 


Ursa Major I 


~100 


~20 


~-2, -2.2''), -2.1^' 


~7'.8 


~-6.5 


Ursa Major II 


32±5 


17.5±0.3 


~-1.3, -1.6'*), -2.0-''' 


~12' 


-3.8±0.6 


Willman 1 


38±7 


17.9±0.4 


~-1.3, -i.e'*' 


~1'.9 


~-2.5 



Note. — Val ues are from the res pect ive discovery papers ( see text), unless other wis e indicated: a ) DaU'Or a et al.l (120061) ; b) 
ISieeell (1200^) : cl IMufioz eFaLl 1 12003) : d) IMartin eFalT (I2007al) : e) llbata et al] l l200l) : f)[S imon Geh a (2007) 



4.1.2. Star formation histories 

To constrain the SFHs of Boo I, CVn I, Leo T and UMa 
II we follow the same approach as before for the GCs. A 
grid composed of nine age bins and twelve metallicity 
bins is used, the age bins going from 10.2 to 8.4 log- years 
with an 0.2 log-year width, and the metallicities from 
[Fc/H]=-2.4 to 0.0 with a 0.2 dex width. A fit is done 
at three different distance moduli, 0.1 magnitude apart. 



with the middle value equal to the literature value. The 
resulting SFHs and AMRs are plotted in Figure [THl 

4.1.3. Object detection 

As with the GCs, we also test our overdensity detection 
scheme, described in Section I2.3.4[ on these new Milky 
Way satellites. The 5°x5° field centered on each ob- 
ject is divided in a square grid of 225 20'x20' boxes, 
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Fig. 13. — Population properties of newly discovered Milky Way satellites: the panels show single component fitting results for Boo 
I, Boo II, CVn I, CVn II and Com, with the columns showing (from left to right): distance vs. metallicity; distance vs. age; age vs. 
metallicity; and the color-magnitude diagram with the best-fit single component model in contours. The gray-scale contours show the 
regions in parameter space that are within Icr, within 2(t, within Scr and outside 3cr of the best-fitting single component model, going from 
white to dark gray. The solid line contours are for the case in which one of the parameters is fixed to the literature value (i.e. age for the 
left-most, metallicity for the middle-left, and distance for the middle-right columns), where available. 



the CMDs of which are fit both with a control field and 
with the control field plus a single component popula- 
tion model. As the control field the whole 5°x5° field is 
used. For the SC model we consider two cases. First we 
use a typical old population model, with age range 10.0 
to 10.1 log-year and [Fe/H]= 2.0 at the best-fit distance 
obtained with the SC fits from Section 14.1.11 Second 
we use a model corresponding to the best-fit properties 
obtained for each object. The former case serves to sim- 



ulate a realistic search of all SDSS data for unknown 
stellar overdensities. When performing such a search, an 
old, metal-poor model would be a likely choice since most 
MW satellites are dominated by such populations. An- 
other factor that would influence the results of a search 
for unknown stellar overdensities is the size of the box 
used. As a test, we also did the same exercise using differ- 
ent box sizes for Boo I. In addition to the 20'x20' boxes, 
we also used 40'x40' boxes and l°xl°boxes. To en- 
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TABLE 4 

Detection significance for new Milky 
Way satellites 





WiQ diaJ I tths 


r ine-tuned uC^ / vttls 


Boo I") 


106.9 


121.6, 58.2, 27.8 


Boo II 


0.86 


0.93 


CVn I 


447.7 


589.4 


CVn II 


2.61 


13.3 


Com 


57.5 


53.9 


Her 


15.2 


28.2 


Leo IV 


13.8 


38.2 


Leo T 


5.8 


11.3 


Seg 1 


9.9 


7.8 


UMa I 


12.9 


11.2 


UMa II 


1156.2 


865.3 


Wil 1 


48.8 


47.8 



Note. — a) the three values of dQ/rms for 
the fine-tuned model fits to Boo I correspond to 
the box sizes of 20'x20' , 40'x40' and l°xl°, 
respectively. 

able the use of these coarser grids, we retrieved a larger, 
10° X 10° field from the SDSS DR5 archive centered on 
Boo I. 

In Figure [T7| we present the results of the fits with 
the fine-tuned stellar population models to the twelve 
objects. Each panel shows the improvement in the 
goodness-of-fit, dQ, as function of the position in the 
field centered on a satellite. Figure [T51 shows the result- 
ing dQ for the three different box sizes in a 10° x 10° field 
centered on Boo I. The significances of all detections, ex- 
pressed in dQ/rms, are listed in Table 21 

4.2. Discussion 
4.2.1. Satellite properties 

We have performed a systematic and uniform numer- 
ical CMD analysis of the recently discovered MW satel- 
lites. How well the properties of an individual object can 
be determined depends on the number of member stars 
present in the CMD and the observed CMD features. 
For some objects, this type of analysis is hard-pressed to 
provide strong statistical constraints due to the sparse- 
ness of their CMDs. In Table [5] we summarize the re- 
sults of our single component and star formation history 
fits. The distance moduli, metallicities and ages listed 
are the unweighted averages of all SC fits that lie within 
1 (7 of the best SC fit and are therefore considered to be 
acceptable fits. The errors are the standard deviations 
of the parameter values of these fits, combined with an 
added systematic uncertainty of 0.1 dex in [Fe/H] and 
0.1 mag in distance modulus to account for isochrone 
uncertainties. In cases where these properties cannot be 
sufficiently constrained we list the metallicities and ages 
we get when assuming the distance modulus from the 
literature, as this is usually the best-determined prop- 
erty. The star formation history results, where available, 
are summarized in Table [5] by the main star formation 
episodes. Below follows an object-by-object discussion of 
these and the star formation history results. 

Boo I. — This dwarf spheroidal is relatively nearby at 
--60 kpc and an RGB, HB and MSTO are aU present 
in the SDSS data, albeit sparsely sampled. The overall 
properties can therefore be determined with relatively 
good accuracy, and the star formation history is the 



best-constrained of all. Especially the distance of Boo 
I is very well constrained thanks to the HB which is de- 
tected at high S/N due to the absence of field stars at 
5 - r <0. The MSTO and slope of the RGB provide 
constraints on the age and metallicity, respectively. Our 
distance, metallicity and age measurements agree very 
well with the findings of other authors based on spec- 
troscopy and deeper imaging data. The star formation 
history shows one significant epoch of star formation in 
the oldest (> 10 Gyr) bin, with no other bins containing 
a significant SFR. Considering the ability of these fits to 
distinguish between a single and a multiple component 
system, demonstrated by the fits to the artificial CMDs 
in Section [2111 we conclude, despite the sparseness of the 
CMD, that Boo I very probably only contains old stars. 
Thus, our results are consistent with the picture of a 
purely old, metal-poor stellar population. 

Boo II. — The number of Boo II member stars in the 
SDSS data is very small and traces only the RGB. Un- 
fortunately, the resulting very low S/N combined with 
the degeneracies between the different parameters make 
it impossible for MATCH to constrain the stellar popu- 
lation parameters in a sta tistically significant way. Al- 
though jwihiEan (|2007l ) show that the SDSS CMD of 
Boo II is consistent with Boo II having similar metallic- 
ity, age, and distance as Boo I, The SDSS data alone do 
not support any significant constraints on any of these 
properties. The only way to get such constraints is by 
obtaining deeper follow-up data. 

CVn I. — Being the brightest of the twelve objects an- 
alyzed here, CVn I has a well-populated RGB and HB. 
This helps to increase the robustness of the constraints 
we obtain from the SC fits. Although at the edge of the 
sensitivity limits of SDSS, the HB enables the determi- 
nation of the distance. Moreover, the extent of the HB 
to very blue colors requires the stars in CVn I to be both 
old and metal-poor. This constraint on the metallicity 
is further enhanced by the slope of the RGB. Together, 
the extended HB and RGB thus can provide robust con- 
straints on the properties of the dominant stellar popu- 
lation. The distance modulus we find is slightly smaller 
than that found in the literature, but not significantly so. 
The metallicity of [Fe/H]=-1.8±0 2 falls niceljy withi n the 
range of metallicities observed bv llbata et al.l (|2006f ). Al- 
though an age of 14±2 Gyr fits the CMD best, the star 
formation history shows two peaks. The oldest star form- 
ing episode in the 10-16 Gyr bin is dominant in terms of 
the number of stars and drives the SC fits, but the peak in 
the 2.5-4 Gyr bin is equally significant. For the old stars 
the metallicity is found to be close to [Fe/H] =-2, while 
the younger stars have metallicities around [Fe/H] =-1.5. 
It should be noted that the young population found here 
is not caused by the algorithm attempting to fit BHB 
stars, as is the case for M15 in Figure |7l Stars belonging 
to this young population only appear in the regions of 
the CMD where the RGB and the red part of the HB 
are located. The presence of two populations is therefore 
purely based on the RGB morphology and the density of 
stars along the HB. This detection of an old and interme- 
diate age popu lation w' i th the se m etallicities confirms th e 
findings of Ib ata et all (l2006l )_and lMartin et all ()2007bf l. 
Furthermore, iZucker et al.l ([2006a) note the presence of 
a carbon star very close to the center of CVn I and its 
heliocentric radial velocity of 36±20 km/s is consistent 
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Fig. 17. — Using MATCH to detect or verify stellar overdensities: Each panel shows (for one of the Milky Way satellites) the difference 
in fit quality, dQ, between fitting only a control field and fitting a control field plus SC model to the CMDs of a grid of boxes on the sky. 
For this, the 5°x5° fields centered on each target were divided in 225 20'x20' boxes. A high value of dQ indicates a large improvement 
and the presence of a large number of stars that are well-fit by the SC model. 



with the radia l veloc ities that llbata et all (|2006l ) and 
iSimon fc Cxehal (|2(¥)1 measure for CVn I, making it is 
likely that this star belongs to the CVn I dwarf. Since 
carbon stars are often of intermediate age, this is an- 
other indication of the presence of an intermediate age 
population in this system. 

CVn II. — Although very faint and distant, our 
distance measurement for CVn II is very accurate 
and in very good a greement with the value from 
iBelokurov et al.l ()2007f ) due to the presence of HB stars. 
From the extension of the HE to very blue colors [g — r^ 
—0.3) we can infer that the stellar population is both old 
and metal-poor. Despite the sparseness of the CMD we 
thus get relatively good constraints on these parameters 
and we find a metallicity of [Fe/H]'^-2 and a very high 



age. CVn II thus appears to be a very old and metal-poor 
object. 

Com. — An absence of a clear HB makes constrain- 
ing the distance of Com difficult and leads to a strong 
degeneracy between age, distance and metallicity. Sim- 
ilar degeneracies occur for artificial CMD 1 in Section 
[lHand the GC NGC 6229 (see Figures [2 and [T]) , even 
though their CMDs are much less sparsely populated. 
The RGB and MSTO alone do not allow a robust deter- 
mination of the individual population parameters. As a 
consequence, the distance modulus we recover has rather 
large error bars but is consistent with the findings of 
IBelokurov eFan (|2007l) . The SC fits favor an ancient 
population (> 10 Gyr) and a metallicity of [Fe/H]~ —2, 
but also here the uncertainties are large. 
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Fig. 18. — Overdensity detection applied to Boo I using different spatial resolutions. A 10° xlO° field centered on Boo I is divided into 
20'x20' (left panel), 40'x40' (middle panel), and l°xl° (right panel) boxes. The fit improvement, dQ, between a pure control field fit 
and a control field plus SC model fit, is plotted as function of position on the sky. 



Her. — In the case of Her the contamination of the 
CMD with fore- and background stars is severe, which is 
due to its location at lower Galactic latitudes compared 
to the other objects. This high contamination means a 
low S/N of the CMD features of Her, apart from the 
HB, which is in the relatively clean part of the CMD at 
g — r < Q. Our distance modulus of 20.4±0.2, which 
largely is dete rmined by the HB is so mewhat lower than 
estimated by iBelokurov et al.l ()2007[ ). but is consistent 
within the error bars. This object too is consistent with 
an ancient and metal-poor ([Fe/H]'~-2) stellar popula- 
tion. Even though the region of parameter space that 
lies within 1 u of the best fit values is well constrained, 
the complete parameter space lies within 2 cr, which is 
due to the low S/N of the CMD. These results based 
only on SDSS data are also consistent with the results 
obtained by applying the same methods to much deeper 
B, V photometry fro m the Large Binocular Telescope 
([Coleman et al.ll2007t ). In fact, the uncertainties on the 
values of distance, metallicity and age from the deeper 
data are very similar to the Icr values we obtain here. 

Leo IV. — The CMD of Leo IV is so sparse that a 
statistically significant determination of the population 
properties is not possible based on the SC fits alone. The 
best fit is given by a population of ~14 Gyr old stars 
with [Fe/H]=— 2.2 at a distance modulus of 20.8, but a 
large number of models gives a goodness-of-fit that is 
only marginally worse. Assu ming a distance modulus of 
21.0 (|Belokurov et al.ll2007l ). gives the solid contours in 
the age-metallicity panel for Leo IV in Figure [TH Using 
this prior, the SC fits allow the age and metallicity to 
be constrained to 14±2 Gyr and [Fe/H]=: -2.1 ±0.3, 
respectively. 

Leo T. — This object is very interesting and differs from 
the others analyzed here, in that it is the least luminous 
system known which shows signs of on-going star forma- 
tion and has HI gas associated with it (Irwin et al. 2007). 
Because of the presence of two clearly very different pop- 
ulations the SC fitting method is not very well suited. 
Even when fitting the two populations separately, by di- 
viding the CMD in two distinct parts as we have done, 
the results are poor, due to the small numbers of stars 
in each part. Taking the results of the two populations 
together we find a distance modulus of 23.1±0.6 mag- 
nitudes. Both the old and young populations are best 
fit with low metallicities, [Fe/H]<-1.5. The star forma- 
tion history we find (Figure fT6|) clearly shows two peaks, 
confirming the presence of two distinct populations. The 
age of the old population, not well-constrained with the 



SC fits, is put in the 6.5-10 Gyr age bin. Comparing 
this with the SFHs of the GCs in Figure [7] and the other 
SFHs in Figurefini where in all cases the peak of star for- 
mation is in the oldest age bin, this implies that the old 
population in Leo T is on average younger than the old 
population in these other objects. The young stars are 
clearly found in the three youngest age bins, where they 
are fit as young (< 1 Gyr) blue loop stars. Ilrwin et all 
(|2007| ) interpreted the nature of these stars similarly. 

Segl.— In Seg 1 only a MS and MSTO are populated 
in the very sparse CMD. As with the artificial CMD 1 
in Section [Ml and the GCs M13, M15 and M92, where 
the SC fits also solely rehe on the MS and MSTO, strong 
degeneracies between age, distance and metallicity are 
expected. Indeed, the contour plots for Seg 1 in Figure 
1141 show a similar morphology, but due to the sparseness 
of the CMD with even poorer constraints. This makes 
any determination of the population properties of Seg 
1, based solely on the SC fits, i mpossible. Even when 
using the distance modulus from IBelokurov et al.l ([2007f ) 
as a prior, the age-metallicity degeneracy prevents strong 
constraints on these parameters, resulting in very large 
error bars on the values for Seg 1 listed in Table [5l 

UMa L — This object also has a very sparsely popu- 
lated CMD. An RGB is readily seen, with some prob- 
able HB stars that, however, are too few to produce 
strong statistical constraints on the distance. Because 
of this and the usual degeneracies, the significance of the 
SC fit results is very low, causing all parameters to be 
poorly constrained. A ssuming a distance modulus of 20 
f Will man etlHI [20051 . our results are consistent with 
an old (^10 Gyr) and metal-poor population, but with 
large error bars. 

UMa IL— Although the MSTO of UMa II is well- 
populated, the scarcity of stars in other parts of the 
CMD provides problems for this object. There are some 
RC stars, but the S/N of this CMD feature is not high 
enough to make for a reliable distance determination. 
Added to that there are strong degeneracies between the 
different solutions; similar degeneracies are seen in Sec- 
tion[2|for the GC NGC 6229 where also the MSTO drives 
the SC fit results. Our results are, however, consistent 
with previous findings. With a prior on the distance 
m odulus, m — M =17.5, the fits agree with the conclusion 
of IZucker et al.l ([2006bD that UMa II is not very metal- 
poor, but seems to have a more intermediate metallic- 
ity of [Fe/H]^ —1.5, although the error bars are large 
(see Table [51). The spe c troscop ic metallicity measure- 
ment from iMartin et al.l ([2007a[ ) is consistent with our 
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TABLE 5 
New MW satellite results 



— — 

Object 




m-M 


[Fe/HJ 


Age 


SF episodes 


Bootes I (Bool) 


18.8±0.2 


-2.2±0.2 


14±2 




-10-16 Gyr 


Bootes II (BooII) 












Canes Venatici I (CVnl) 


21.5±0.2 


-1.8±0.2 


14±2 


-2-4 Gyr, r 


-10-16 Gyr 


Canes Venatici II (CVnll) 


20.8±0.2 


-2.1±0.3 


14±2 






Coma Berenices (Com) 


18.4±0.4 


-1.9±0.4 


11±5 






Hercules (Her) 


20.4±0.2 


-2.1±0.2 


14±3 






Leo IV (LeoIV) 






14±2 






Leo T (LeoT) 


23.1±0.6 


-1.6±0.6 




< 1 Gyr, 


-6-10 Gyr 


Segue 1 (Segl) 




-1.6±0.5 


n±4 






Ursa Major I (UMal) 




-1.8±0.4 


10±5 






Ursa Major II (UMall) 




-1.5±0.5 


ll±4 


~2-6 Gyr, - 


-10-16 Gyr 


Willman 1 (WiH) 


17.9±0.4 


-1.6±0.5 


10±5 







value, but the one from iSimon fc Gehal (|2007h IS more 
metal-poor, but still consistent within the errors. While 
the SC fits imply that the majority of stars is likely old, 
~11 Gyr, the SFH (Figure [T6|) hints at more complexity. 
According to the SFH results, there is a population of 
very old (>10 Gyr) and metal-poor ([Fe/H]~-2), as well 
as a younger ('^5 Gyr), close to solar metallicity popu- 
lation. Since the SC fits measure the average properties, 
this is con sistent with the inter mediate metallicity we 
find there. IZucker et al.1 (|2006b[ ) report that the upper 
main sequence of UMa II in their deeper Subaru data is 
wider than expected just from observational errors, one 
explanation for which could be the presence of a mix of 
populations. Also in our fit, the double peak results from 
the width of the MSTO, which is less well-fit by a single- 
age population. The exact ages and metallicities of the 
populations should be studied with deeper and spectro- 
scopic data, since the SDSS CMD does not support such 
detailed conclusions. However, taking our SFH fit and 
these other results together it does seem likely that UMa 
II, like CVn I, is a system with a complex star formation 
history, which might also be (part of) the explanation for 
the discrepancy between the two available spectroscopic 
metallicity measurements. 

Will. — For this system only the MSTO is observed, 
which is also sparsely populated. This leads to a strongly 
degenerate solution space, similar to the cases of UMa 
II and Seg 1. The distance can be determined slightly 
better here, but the range of distance moduli for which 
fits within la from the best fit can be obtained is still 
very large: 1 7. 2< m — M <18.5. U sing the distance 
modulus from IWillman et all (|2005bD (m-M=17.9) as a 
prior, the fit results imply an old population ('^10 Gyr) 
with an intermediate to low metallicity ([Fe/H]'^-1.6), 
but also here the error bars on these numbers are large. 

Figure [H] summarizes our results from the MATCH 
fits to the CMDs for the sample of new Milky Way satel- 
lites. In the left panel we plot our metallicity estimates 
for all objects versus their absolute magnit udes, together 
with th e sp ectroscopic met a llicitie s from iMartin et al.l 
(|2007aD and ISimon &: Gehal (|2007[ ). Especially for the 
fainter objects our metallicity estimates are very poorly 
constrained (Figures [T3] and [T3]) and because of the avail- 
able isochrones our metallicity range is restricted to 
[Fe/H]= —2.4 at the metal-poor end. In light of this, 
our photometric results agree remarkably well with the 
spectroscopic values. Still, care has to be taken when 



drawing conclusions from these metallicity estimates. We 
can say at best that we find no indication for a corella- 
tion between metallicity and the luminosity of the ob- 
jects, which is observed for bright e r dwarf spheroidal s 
(e.g.. lGrebel & Guhathak^flOM iGrebel et al.ll200l . 
The mean age of the stars versus absolute magnitude is 
shown in the middle panel of Figure 1191 For objects in 
which the SFH fits showed evidence for two populations, 
both populations are plotted separately. The right panel 
shows the inferred ages of the satellites versus distance 
modulus, to probe possible environmental infiuences. 

Looking at the middle panel of Figure [19] there is no 
clear evidence for a luminosity dependence of the abil- 
ity of a dwarf galaxy to form stars recently. Although 
the two most luminous objects in this sample, Leo T and 
CVn I, both have young or intermediate age populations, 
the much fainter UMa II does as well. While the environ- 
ment might play a role, the right panel of Figure [T9l shows 
no obvious relation between age or the presence of mul- 
tiple star formation epochs with distance from the sun 
(since these objects are at high Galactic latitude and at 
large distances, the heliocentric and galactocentric dis- 
tances are not significantl y different). T h is is contrary 
to the trend suggested by Ivan den Berghl ()1994( ) for the 
brighter Local Group dwarf spheroidals, in which sys- 
tems with a larger fraction of young or intermediate age 
stars are systematically more distant. 

4.2.2. Satellite detection 

Since most of these recently detected objects are in- 
trinsically faint and/or distant, they are of course much 
harder to detect by any overdensity detection scheme 
than the GCs. The improvement in the goodness-of-fit, 
dQ, between a fit with a control field only and a con- 
trol field plus SC model is roughly proportional to the 
number of stars in the CMD that can be fit by the SC 
model. This means that the objects with sparsely popu- 
lated CMDs are the most difficult to detect. Comparing 
the overdensity detection results in Figure [T7] with the 
CMDs of the objects (Figure [12]) confirms this, as Boo, 
CVn I and UMa II give by far the highest dQ. Similarly, 
the rms noise in dQ increases with increasing numbers 
of field stars in the CMDs, and at the low S/N end the 
noise level starts to play an important role. For example, 
comparing the panels for Boo II and CVn II in Figure 
[T71 which have the same dQ scale, reveals very different 
noise levels. This can be traced back to the positions of 
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Fig. 19. — CMD fitting results for new Milky Way satellites: Left: metallicities from this paper (filled ci rcle s) versus absolute magnitude, 
compared with spectroscopic metallicities fromjMartin et al. (200f3) (open squares) and Simon & Geha ( 20075 (open triangles). For clarity, 
only the error bars on the results from this work arc shown. Middle: ages of the satellites versus absolute magnitude. Open circles denote 
objects for which a full SFH fit was done, and in the cases where two distinct populations were found (CVn I, Leo T and UMa II) both are 
plotted; filled circles denote objects where the age is only determined from the SC fits. Right: as in the middle panel, the ages are plotted, 
but now versus distance modulus. 
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Fig. 20. — Histogram of the dQ/rms distribution for an empty 
region of sky of effective area of ^SOO square degrees. 



these objects on the sky, with Bootes II lying behind the 
Sagittarius stream and CVn II in a much 'cleaner' region 
of the halo. The Sagittarius stream also produces a high 
noise level around Segue 1. 

While most objects produce clear peaks in dQ, the de- 
tectability of several is rather low. For example, CVn II, 
Leo T, Segue 1 and UMa I have a dQ only a factor of a 
few higher than the most prominent noise peaks; Boo II 
is not clearly detected at all. When inspecting the values 
in Table |4l it turns out that the use of a fine-tuned popu- 
lation model as compared to the generic old, metal-poor 
model has a varied effect. In some cases the detection 
significance is improved, in some cases the change has 
little effect and in a few cases the detectability is actu- 
ally decreased. In these last cases, even though the peak 
dQ is the same to or higher for the fine-tuned model, the 
change in the rms causes the detection to become less 
significant. This indicates that, while sampling a range 
of population models would be preferable when doing an 
overdensity search, a relatively small grid of models will 
not severely compromise the outcome. 

To test the significance of the dQ/rms measure, we 
tested the overdensity detection scheme on a presum- 
ably 'empty' patch of SDSS sky with 170°<RA<200°and 
40°<Dec<65°. This patch of sky was divided into 
20'x20' boxes with some overlap, giving a total of 7088 
boxes, equalling an effective search area of almost 800 



square degrees. A control field (the CMD from the 
surounding 5°x5°area) was fit to the CMD of each box, 
as well as this control field plus a typical old stellar popu- 
lation model with an age bin of 10.0 to 10.1 log-year and 
[Fe/H]==-2. This was done seven times, with the stellar 
population model shifted to different distance moduli: 
17.0, 18.0, 19.0, 20.0, 21.0, 22.0 and 23.0 magnitudes. 
Thus, for each box seven values of dQ/rms are obtained, 
where the rms is obtained from the boxes in the suround- 
ing 5'^x5°. The methodology closely resembles the ap- 
proach used for the MW satellites. Figure [20l shows the 
histogram of all dQ/rms values for all boxes (solid line), 
as well as for two individual distance moduli, 18.0 (dot- 
ted) and 22.0 (dashed) mag. All three histograms show 
the same distribution, implying that the noise behaves 
similarly for all distance moduli. The tail of the distri- 
bution extends to dQ/rms ^ 8. An inspection of the 
locations with values higher than 7 reveals no obvious 
overdensities, implying that this tail is indeed caused 
by noise. In order to obtain a clean sample of stellar 
overdensity detections, a dQ/rms threshold of ^^9 seems 
appropriate. This means that when using a sufficiently 
sampled grid of population models and distances, all new 
MW satellites would be detected, with the exception of 
Boo II (cf. Table HI). 

The results of the fits to Boo I with different box sizes, 
listed in Table H] and shown in Figure [HI illustrate the 
effect of the spatial resolution of the CMD grid. Here, 
the detectability is also strongly affected by the rms of 
the background noise. Since the dQ due to random vari- 
ations in the field population increases with increasing 
numbers of field stars, the rms is higher for larger box 
sizes, favoring the smallest boxes used here. The fact 
that the AQ/rms decreases roughly by a factor of 2 as 
we move to larger box sizes is caused by the number of 
background stars increasing by a factor 4, while the num- 
ber of Boo I stars in the central box does not change. It is 
also interesting to note that here the noise is dominated 
by the Sagittarius stream, which runs accross the lower 
left corner of the panels in Figure (TH] The significance of 
the detections in Table |4] can be boosted by decreasing 
the box sizes, especially for the objects with the smallest 
(projected) size. 
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5. SUMMARY AND CONCLUSIONS 

We have shown that apphcation of MATCH (jPolphinl 
12001') to Hess dia grams or CMDs of SDSS data success- 
fully recovers the known stellar population properties of 
globular clusters (Section[3]), within the limitations of the 
data. This gives confidence that despite inaccuracies in 
the theoretical isochrones in isolated parts of the CMD, 
the overall CMD-fitting technique still provides good and 
trustworthy results. Apart from the standard star forma- 
tion history fitting method, we have also demonstrated 
three other ways of applying MATCH to SDSS data: 

1) Identifying the best fitting 'single component' stellar 
population, a mode useful for CMDs of very faint or dis- 
tant objects with limited S/N in their CMD; 

2) Fitting the line-of-sight distribution of stars - to keep 
the fit's degrees of freedom constant, we fix the metallic- 
ity at each age bin; and 

3) Algorithmic detection of objects, i.e. of a localized 
stellar population overdensity in a particular (unknown) 
direction and at a particular (unknown) distance. This 
is done by comparing the difference in the goodness-of-fit 
between a CMD-fit to a control field and a CMD-fit to a 
control field plus model population. 

These tests show that the application of CMD-fitting 
techniques to wide-field survey data is a very promising 
endeavor in the light of SDSS and other surveys that 
will start yielding data over the next few years. Further 
improvements in software and theoretical isochrones can 
only increase the use and accuracy of such work. On the 
basis of these tests and the application to known globu- 
lar clusters, we have applied our single component fitting 
technique to twelve recently discovered Milky Way satel- 
lites. The sparseness of several of the CMDs of these 
objects inhibits a robust determination of the popula- 
tion parameters and the best fit solutions are often not 
unique. In order to obtain good constraints at least some 
CMD features such as the MSTO, HB and RGB should 
be observed with sufficiently high S/N. For the least lu- 
minous objects this is inherently problematic because of 
their lack of (evolved) stars. For only four of the satel- 
lites, the CMDs had sufficiently high S/N to warrant a 
star formation history fit. 

The results of the population fits to the new Milky Way 
satellites, summarized in Tabic [SJ are consistent with the 
individual discovery and follow-up studies. When look- 
ing at the sample as a whole, we find that the majority of 
these objects are likely dominated by old (> 10 Gyr) and 
metal-poor ([Fe/H]~-2) stellar populations. For Boo, 
CVn II, Com, and Her the single component fits converge 
to this picture, and for Boo the full fit of a star formation 
history also confirms this. Seg 1, UMa I, UMa II, and 
Wil 1 arc not well constrained by the SC fits, but using 
the distances from the literature as priors, they are also 
consistent with similar ages and metallicities, although 
slightly higher metallicities ([Fe/H]^ —1.6) cannot be 
excluded. The single component fits show CVn I to be 
old as well, but not extremely metal-poor. However, the 
SFH reveals the presence of both an old, metal-poor, as 
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well as a younger (<3 Gyr) and more metal- rich popula- 
tion. This is consistent with the findings of llbata et al.l 
(|2006[ ). who detect two kinematically distinct popula- 
tions in CVn I which differ in age a nd metallicity. and the 
detection of a young component bv Marti n et al.l ()2007bl ) 
in deeper photometry. The SFH of UMa II also shows 
a peak at an age of ~5 Gyr, which suppor ts suggestion 
of mult iple populations in this object by IZucker et al.l 
(|2006bD . who find the MSTO width of UMa II to be 
too large to be explained by photometric errors alone. 
Finally, in Leo T the presence of two distinct popula- 
tions is readily seen in the CMD and reproduced in the 
star formation history. The RGB seen in this system 
belongs most likely to an intermediate age (~5-10 Gyr) 
population with a metallicity of [Fe/H]~-1.5 to -2, while 
the blue stars in the CMD are most likely very young 
(<1 Gyr) blue loop stars with similar metallicities. The 
objects studied here are all intrinsically faint and sparse, 
and in some cases the population parameters can be con- 
strained only poorly. However, their star formation his- 
tories do not appear uniform, nor do all of them appear 
to be exclusively old and metal poor. 
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